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I. INTRODUCTION 

With the recent discovery of the Higgs boson by the LHC m, the Standard Model (SM) 
has reaffirmed itself as the theoretical explanation of elementary particle dynamics, including 
presenting a plausible picture of Electroweak Symmetry Breaking (EWSB) [5}[B] through the Higgs 
mechanism. While this is so, and not withstanding the fact that SM has been extensively tested 
very successfully by many different experiments [7], it is widely believed that the SM is an effective 
theory, which needs to incorporate suitable modifications at large energies. Many expect that this 
large energy scale could be as small as a few TeV, which is being explored at the LHC. Apart from 
direct measurements of the properties of the Higgs boson, like its interaction couplings with itself, as 
well as with gauge bosons and fermions, signature of EWSB could be probed by understanding the 
structure and values of gauge-boson self interactions. This is so, because the longitudinal degrees 
of freedom of the gauge bosons arise from the Higgs sector. The relation between the Higgs-gauge 
boson couplings and the trilinear gauge couplings are presented and studied in Ref. mm- The 
discovery of Higgs boson of about 126 GeV mass, while establishing the Higgs mechanism as the 
method of EWSB, has opened an era of Higgs precision studies. Being a hadronic machine LHC 
has limitations to undertake precision studies. The International Linear Collider (ILC) proposed 
to collide high energy, high luminosity electrons and positrons has the mission of studying the SM 
at high precision and to look for signals beyond the standard model [10] [Tl]. Such a machine is 
well suited to an in-depth analysis of the gauge boson interactions within and beyond the SM. A 
large number of physics studies establish the fact that the potential of ILC is further enhanced by 
considering high energy photon-photon (77) collisions as well as electron-photon (e7) collisions, 
apart from the electron-positron (e + e _ ) collisions. 

It is obvious that, the e7 and 77 colliders are better suited to study the photon couplings 
with other gauge bosons like the 7WW, yylTW, 7 ZZ and 77 Z [12]. In this article we will 
focus on ey —» vW with W —► Iv. Signature of such an event is a single lepton with large 
missing energy. This process is sensitive to new physics effects including anomalous 7 WW USE], 
composite fermion models mi, etc. In the case of anomalous 7 WW couplings, this process has 
the advantage over e + e _ —► W + W~ [TS], which is sensitive to both 7 WW and WW Z couplings. 
Again, 77 — > W + W~ [T7, jT3] is sensitive to 7 WW, 77 WW and 77 Z couplings, revealing the edge 
of ey collider to study -)WW. In most of the previous studies, observables at the production level 
of the single W is investigated, with the exception of Ref. m, where the authors consider angular 
spectrum of the secondary leptons, also including the effect of off-shell W. Our main aim of this 
work is to study the possibility to exploit the secondary lepton spectrum including the energy and 
angular distributions to probe relevant new physics signals arising in anomalous 7 WW coupling. 

In many models beyond the SM the quartic and triple-gauge boson couplings including 7 WW are 
altered from their SM values. In a model independent approach, an effective Lagrangian with terms 
additional to the SM Lagrangian is considered in phenomenological and experimental studies m- 
Relevant to the process considered here, the effective *yWW vertex is commonly parametrized in 
terms of Sk 7 and A 7 , in the absence of CP violation (with vanishing SM values). LEP constraints in 
single-parameter analysis (taking one parameter at a time, keeping the others at their SM values) 
gives bounds of —0.105 < Sk 7 < +0.069, —0.099 < Sk 7 < +0.066 and —0.059 < A 7 < +0.026, 
—0.059 < A 7 < +0.017 and two-parameter analysis limits their values to —0.072 < i5k 7 < +0.127 
and —0.068 < A 7 < +0.023 [50;, ]5T| at 95% C.L. Tevatron constraints from Wy process are not 
contaminated by other couplings, but are more relaxed compared the LEP constraints to give 
—0.51 < 6 k 7 < +0.51 and —0.12 < A 7 < +0.13 [22]. The CMS constraints [23j from Wy process 
are also more relaxed compared the LEP constraints. 

Phenomenology of anomalous ^WW coupling in the context of LHC as well as ILC has been 
carried out in a number of recent publications [13 E3 WMM- In particular, m has analyzed single 
W production with its leptonic decay to probe the effect of anomalous couplings in ey collision. 
In this work we exploit the full potential of the secondary lepton spectrum to study the effect of 
7 WW coupling in ey -+ vW —> v(lv). 

In the next section we provide some details of the process and the observables used. In section 
3 we present our numerical results, and finally summarize the study and present our conclusions 
in the last section. 
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II. ANALYSIS AND DISCUSSION 


Considering a real on-shell photon, the most general CP-conserving 7 WW coupling within a 
Lorentz invariant Lagrangian can be written in the following form m- 


c 


7 WW 


W^W^A" - + (1 + <k 7 ) W$W v F> tv + W^W? F x 


°W 


(i) 


In the SM, the gauge structure SU(2)l x U(l)y dictates the jWW couplings, with vanishing 6k 7 
and A 7 at tree level. Therefore, precise measurements of these couplings will test the gauge sector 
of the electroweak interactions. Fig. 0 shows the Feynman diagrams for the process along with the 
momentum labels. 



FIG. 1: Feynman diagrams for the process ey —¥ vW —> First diagram shows the momenta used 

and second diagram shows the indices used in the vertex in Eq. [3] 


With the effective Lagrangian given by Eq. 0 the 7 WW vertex for the process under study 

(Fig-0 (b)) takes the form: 


f^a = ie 


2A4 g" x + 2k" g» x - (fc 7 + k w ) x g» v + (fe 7 - X^k" g» x - k x g» v ) (2) 


+ + k w ) x (k^k" - (fc 7 • k w )g^) 

m z w 


The effective Lagrangian in Eq. 0 should be considered as a low energy approximation of some 
fundamental theory, which is expected to emerge at some high energy scale, A. To control unitarity 
violation at high energies, we consider the anomalous couplings as form factors according to the 
following M 


A = A 0 / 




|(fc 7 - k w ) 2 \\ 
A 2 ) 


(3) 


where Aq = 5« 7 , A 7 . 

Contribution from the third diagram (Fi S0> is found to be a few percent, and will be neglected. 
This allows us to perform our computation in the Narrow-Width Approximation (NWA) in which 
the W-propagator is approximated to get 


1 

\ k W — m w \ 2 


7r 

mw^w 


${ k w ~ m w\ 


(4) 


where mw is the mass and Tw is the width of the W boson. 
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FIG. 2: Reference frame defining different momenta and corresponding angles used in Eq. [5] 


To perform phase space integrations we fix our reference frame as the center-of-mass frame 
(CMF) of the colliding electron and photon system, z-axis is taken along k^, which is the momen¬ 
tum of the outgoing lepton (considered as muon in the further discussion) as pictured in Fig. [2] 
y-axis is defined as k M x k e , where k e is the momentum of the colliding electron. The v M comes 
out at a polar angle and azimuthal angle . Energy-momentum conservation and the NWA 
(Eq. [4]) are used to get the differential cross-section 


da 

dx^ dcos 0 M d<pj? fl 


(27r) 3 128 rriw r w 


| M r 


( 5 ) 


where x M = -j*, 


with y/S the center-of-mass energy and E^ is the energy of the muon, and 


cos 9^ 
as, 


kjj.-ke 

I -fee | 


Here M r is the reduced amplitude given in terms of the invariant amplitude M 


M = 



M r . 


( 6 ) 


\M r \ 2 is obtained using FORM computational package [27] . After integrating the unobservable 
, we get the double distribution of energy and polar angle of the secondary muon in the center- 
of-mass frame of the colliding particles with the electron momentum now taken along the redefined 
z axis. Notice that the muon energy in the CMF is bounded by < E M ^ ■ 

To obtain the distribution in the lab frame, we need to boost the above differential cross-section 
appropriately. For an electron beam energy of E e and photon energy of ui 1 = x E e , we have the 
following relation between the variables in the CMF and the laboratory frame. 

s = x 4 E 2 


E, 

cos 9 M 


K ab 7 (I-/? cos# 6 ) 

cos # b - ft 
1 — /? cos ’ 


( 7 ) 


where /3 is the speed of the CMF compared to the lab frame, and 7 = 




limits of £3 ab integration depends on cos # , keeping it within the bound 


Notice that the 


m w _ 1 < E iab < Vx E e 

4y/xE e 7 (l — fd cos # 6 ) — ^ 7 (l — /3 cos # b ) 


( 8 ) 
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It is well known that the hard photons produced by compton scattering are polarized. The po¬ 
larization depends on the initial electron and laser beam polarizations. Besides, the final photon 
distribution will depend on the initial beam energies. Considering these, the parton level cross sec¬ 
tion need to be folded with the appropriate luminosity function, C ei (x ), an expression for which 
is provided in the Appendix, so that the total cross-section in the lab frame is given by 


a = 


d/L*yj e(*T-) 
dx 


a(s) dx. 


(9) 


We suitably use Eqs. 00 EH to obtain the total cross-section, the muon angular distribution 
rfcose iai > > the scaled muon energy distribution c // ab , and the energy-angle distribution of the muons 

dx ,ab d cos o lab the t & h f rame - Integrations over the photon distribution variable x, the scaled muon 

energy x l / b and the muon angular variable cos 9 l / b in appropriate cases are performed numerically 
using the Cuhre routines under the CUBA package [28j . From now on we will drop the superscript 
lab from the variables. 

Phenomenological analysis of ey —>• vW is considered by some authors in the past mm- 
Most of the studies limit their analysis at the production level. Experimentally it is more useful 
to understand the effect on the final state particles arising from W decay. This is all the more 
important in the case of leptonic decay, as it is not possible to reconstruct the W in such case. 
In ref. |14) . analysis including decay spectrum is presented, where detailed study of the secondary 
lepton angular distribution is considered along with other reconstructed observables concerning 
the W production. In the analysis presented here we demonstrate the usefulness of the combined 
energy-angle distribution of the secondary leptons in extracting information on the anomalous 
7 WW couplings. While Ref. m takes into account the contribution due to off-shell W along 
with the on-shell production, our work is in the NWA assuming on-shell W production. At the 
same time, as note above, the effect of off-shell contribution, mainly arising through the diagram 
with muon propagator, is negligible, while the our analysis with energy-angle double distributions, 
along with combinations of the angular and energy distributions gives a handle on disentangling 
the effect of c)re 7 and A 7 , to certain extend. 


III. NUMERICAL RESULTS 

For our numerical analysis we consider ILC of center-of-mass energies of 500 GeV and 1000 GeV, 
with the option of ey collision using backscattered laser photons, and with unpolarized electron 
beams. It may be noted that as the process being considered is a purely a weak interaction process, 
only left-handed electrons will take part. Therefore, the left-polarized electron beam would enhance 
the cross-section, without presenting any additional advantage. We shall present all our analysis 
for unpolarized beam, without losing any generality. 

Firstly, we consider the total cross-section. Fig. [3j shows the total cross-section against the 
center-of-mass energy of the e _ e _ system (denoted as y/s), one of which Compton-scatters on the 
laser beam to produce the high energy photon beam. The effect of <5 k 7 has an effect of about 7% 
to 8 % almost all through the different center-of-mass energy values, while the effect of A 7 has a 
dependence on the y/s. This is expected, as the former is the coefficient of a dimension-4 operator, 
while the latter is the coefficient of a dimension -6 operator, thus producing a momentum depen¬ 
dence in the corresponding coupling. Again, for the same order of magnitudes of the couplings, the 
cross section is more sensitive to <5 k 7 compared to A 7 , as expected for higher dimensional opera¬ 
tors. One may note that the sensitivity to the sign of A 7 ceases beyond the center-of-mass energies 
slightly beyond 500 GeV. This means that the effect of the interference term is overshadowed by 
the contribution purely coming from the dimension -6 anomalous term, which goes like the square 
of A 7 . In order to be more quantitative, in Table |T[ we have presented the cross-sections for the 
limiting values of the couplings for two selected center-of-mass energy values, y/s = 500 GeV and 
1000 GeV. Taking individually separately, the effect of A 7 is small, with a maximum of about 7% 
for negative values at 1000 GeV center-of-mass energy. It may be noted that the effect is sensitive 
to the sign of the parameter, indicating that the interference term is dominant, as expected. At 
500 GeV, one can ignore the effects of A 7 even when considering along with <S/c 7 , whereas at 1000 
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Vs [GeV] 


FIG. 3: Total cross-section against yfs = 2 E e , where E e is the electron beam energy, for different 
anomalous couplings, along with the SM case. 




FIG. 4: The surface plots showing the total cross-section against 8k 1 and A 7 for the center-of-mass energy 
(y^s) of 500 GeV (left) and 1000 GeV (right), where cross-section in pb is along the vertical axis. 


GeV, the presence of A 7 can either nullify or enhance considerably the effect of i5k 7 . Largest effect 
is seen when i5k 7 is positive and A 7 is negative, with about 20% deviation in the total cross-section 
at 1000 GeV. In the surface plots shown in Fig|4j variation of the cross-section with <5 k 7 and A 7 
at center-of-mass energy of 500 GeV and 1000 tieV are presented, which show monotonous de¬ 
pendences, almost a linear one, as expected from the fact that only the interference terms play 
role. We shall see below that, with such large cross-section, small deviations even at percent or 
sub-percent level will be able to probe the relevant parameters significantly. 


The two-parameter bound showing the allowed region in the d« 7 — A 7 plane is presented in Fig[5j 
at center-of-mass energy of 500 GeV with an integrated luminosity of 100 /6 _1 .This corresponds 
to number of event of about 430000. Thus, with the dependence described above, the reach on the 
parameters for these machine parameters is expected to be very good. We remind the reader that, 
when cross-section is considered as a function of anomalous coupling parameters <5 k 7 and A 7 , the 
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5 K>ry , 

A'y 

a [pb] 



yds = 500 GeV 

y/s = 1000 GeV 

SM 

case 

4.284 

4.568 

+0.066, 

0 

4.626 

4.977 

-0.072, 

0 

3.922 

4.135 

0, 

+0.017 

4.319 

4.602 

0, 

-0.059 

4.238 

4.942 

+0.066, 

+0.017 

4.662 

5.012 

+0.066, 

-0.059 

4.576 

5.350 

-0.072, 

+0.017 

3.955 

4.169 

-0.072, 

-0.059 

3.879 

4.510 


TABLE I: Total cross-section for different combinations of TGC parameters at \fs = 500 GeV and 1000 
GeV. 


result is a second order polynomial in these two parameters. With this, the 3er limit of the cross- 
section leads to an elliptic equation corresponding to the relation between these two parameters. 
This results in an elliptic band in the 5k 7 — A 7 plane respecting the 3 a limit of the cross-section. 
Very strong bound on 5k 7 can be derived for any assumed value of A 7 . On the other hand, for 
small values of i5k 7 the sensitivity of cross-section on A 7 is comparatively smaller. At the same 
time, the relative sign of the parameters is preferred to be opposite to each other, as was also 
evident from the analysis above. 



FIG. 5: The shaded region corresponds to values of <5 k 7 — A 7 with the total cross-section within the 3<r 
limit, for an integrated luminosity of 100 fb ~ 1 at a center-of-mass energy of 500 GeV. 


Moving on to the single parameter, sensitivity of cross-section on one parameter, while the 
other one is fixed is presented in Fig. [Gj for center-of-mass energies of 500 GeV and 1000 GeV. 
The green band represents the 3<r region of the SM value of the cross-section, considering an 
integrated luminosity of 100 fb -1 . The cross-section has a linear dependence on <5 ac 7 , showing 
that the contribution proportional to the quadratic term is negligible. The 3er allowed range of 
5k 7 shifts with the value of A 7 . Although, as Table |T] suggests, the effect of A 7 is small compared 
to the large effect of 5k 7 values, the influence is significant for derivable limits as can be read 
from Fig. [6] For example, at 500 GeV center-of-mass energy, assuming A 7 = 0, one obtain a 3 a 
limit of —0.004 < 5k 7 < +0.004, which is moved to about —0.008 < 5k 7 < 0 for A 7 = 0.01, or 
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FIG. 6: Cross-section against 5k 1 (top row) and A 7 (bottom row), when the other parameter assume 
typical values. The center-of-mass energies considered are yfs = 500 GeV (left column) and 1000 GeV 
(right column). The green band indicates the 3a limit of the SM cross-section, with integrated luminosity 
of 100 fb' 1 . 


moved to 0 < 5k 7 < +0.008 for A 7 = —0.014. In the case of 1000 GeV center-of-mass energy, 
the cross-section is not sensitive to the sign of A 7 , as already discussed. Thus, the deviation of 
the dependence of the cross-section of A 7 from linearity is evidently seen in Fig. [6] (lower right). 
On the other hand, in the case of A 7 , at 500 GeV with 100 fb _1 , and keeping dre 7 = 0, a limit of 
—0.013 < A 7 < +0.01 could be reached, which is changed to —0.02 < A 7 < 0.004 if 5 k 7 = +0.002, 
and —0.005 < A 7 < 0.015 if 5 k 7 = —0.002. In the case of 1000 GeV center-of-mass energy, the 
reach on A 7 are —0.013 < A 7 < +0.013, —0.008 < A 7 < +0.008 and —0.017 < A 7 < +0.017 
corresponding to Sk 7 = 0, + 0.002 and —0.002, respectively. Thus, quite evidently, the study of 

the cross-section alone will not provide information regarding the values of or reach of either of 
the parameters in question. Possibilities of disentangling these effects will be discussed through a 
study of different kinematic distributions of the decay products. 

We next consider the angular distribution of the secondary muons plotted in Fig. [7] for different 
combinations of 5 k 7 and A 7 , at center-of-mass energies of 500 GeV and 1000 GeV. As can be seen, 
the effect is more pronounced at smaller cos 9 M values, while most of the events are gathered in the 
forward direction. With high luminosity, one can expect high statistics for single W production 
at ey collider. This means, that either the entire backward region, or even a smaller region of the 
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FIG. 7: The angular distribution for different combinations of 5k-, and A 7 for the two different values of 
y/s = 500 GeV (left) and 1000 GeV (right). 


cos 6^ could be probed to see the effect of anomalous TGC. In Table [TT] the deviation from SM case 
is illustrated by considering a bin of —0.90 < cos 6^ < —0.85 with an integrated luminosity of 100 
fb^ 1 . The number of events in the backward hemisphere (7V back ), and the deviation from the SM 
case (A back = V back — N^ k ) show that deviations as big as 30% is possible in at yfs = 500 GeV, 
which is increased to close to 150% at yfs = 1000 GeV. The deviations are more pronounced in the 
case of 1000 GeV, compared to 500 GeV. The large sensitivity of this observable to the different 
combinations of the two anomalous couplings gives a way to distinguish different scenarios. A very 
moderate luminosity of 100fb _1 gives enough statistics to mean that 10% deviation correspond to 
5 - 6 a deviation. Another way to present this effect is through the forward-backward asymmetry, 
defined as 

f-i rav dcos( ^ " Jo 1 d£e^ dcoB0 i> 

A a — ± _____ ±_ 

— r 1 5 

I-i j^ dcos + fo j^h; dcos 


I \Ano. _ aSM I 

FB (%) = 1 FB a sm FB ' x 10 °- (!0) 

a fb 

In Fig. [8j this asymmetry is presented for different combinations of the parameter values. In the 
SM, at 1000 GeV with 100 fb -1 integrated luminosity, we expect about 460000 number of events in 
total (with 100% efficiency of detection). Out of this, about 380000 is expected to be in the forward 
direction giving rise to an asymmetric number of events of about 300000. Considering the statistical 
uncertainty, even a deviation at a few percent level is perceivable. We have considered a typical 
1% deviation to find the reach of the couplings for this configuration of the collider. The green 
band shows the region that falls within 1% deviation from the SM values. The single parameter 
reach on 5k-, is better through the consideration of cross-section compared to the forward-backward 
asymmetry. On the other hand, the case of A 7 is quite different. At y/s = 1000 GeV, the reach 
could be improved by factor of 2 to —0.005 < A 7 < +0.005. More importantly, the dependence of 
the derivable limits, which was about 30% in the case of the limits drawn from the cross-section, 
is now reduced to about 10% at this center-of-mass energy, clearly demonstrating the advantage 
of the Aps in probing the influence of A 7 . 


So far, we have noted the influence of one of the couplings in deriving the bound on the other. 
In order to see the possibility of disentangling the effects of these couplings, we shall turn to the 
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5k~! 

A-y 

yfs = 

500 GeV 

y/s — 

1000 GeV 



^yback 

A back % 

^yback 

A back % 

SM 

case 

2900 

0 

1600 

0 

0.066 

0 

3250 

12 

2100 

31 

-0.072 

0 

2600 

-10 

1100 

-31 

0 

0.017 

2800 

-3 

1250 

-22 

0 

-0.059 

3400 

17 

2850 

78 

0.066 

0.017 

3100 

7 

1700 

6 

0.066 

-0.059 

3750 

30 

3900 

144 

-0.072 

0.017 

2400 

-17 

750 

-53 

-0.072 

-0.059 

3050 

5 

2350 

43 


TABLE II: The number of events within —0.90 < cos 8^ < —0.85 for different combinations of 5 k; 7 and A 7 
at y/s = 500 GeV and 1000 GeV, along with the corresponding deviation from the SM case. An integrated 
luminosity of 100 fb -1 is considered. 



FIG. 8: The forward-backward asymmetry againstdfv 7 (top row) and A 7 (bottom row) at yfs = 500 GeV 
(left) and 1000 GeV (right). The green band correspond to a ±10% deviation from the SM case. 


kinematic distributions like the energy distribution, alone and in combination with the angular 
distribution. In Fig 9 the energy distribution against where E e is the electron beam 

energy is presented, snowing the effect of anomalous couplings for different combination of TGC 
parameters at yfs of 500 GeV and 1000 GeV. Once again, we emphasis that it is only through a 
combination of different observables that one may be able to disentangle information regarding the 
and A 7 couplings. The influence of Sk 7 remains more or less the same throughout the range of 
x^. At the same time, A 7 > 0 has a slight diminishing effect in the low energy region, and a slight 
enlarging effect at the high energy region. On the other hand, A 7 < 0 has a large enlarging effect 
at low energies , which turns the other way around, but with a much smaller magnitude at higher 
energies. Thus, the energy distribution of the secondary muons become a clear discriminator of 
these qualitatively different scenarios, as summarized in Table |Hl| 


A more detailed picture could be made available by exploiting the energy-angle double distribu¬ 
tion, as presented in Fig. [lOj where the SM distribution (top left) and the deviations from the SM 
in the presence of <5 k 7 and A 7 , individually, are shown. Deviation in the case of A 7 > 0 is negligible, 
and therefore not presented. Evidently, in order to probe A 7 , one need to focus on the low energy 
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FIG. 9: The energy distribution of the secondary muon for different combinations of <5 k 7 and A 7 , 

compared with the SM case. Center-of-mass energies considered are y/s = 500 GeV (left) and 1000 GeV 
(right). 


Scenario 

Effect on ^ energy 

<5k 7 > 0 
A 7 > 0 

enhancing effect 

both low and high energies 

<5k 7 < 0 
A 7 < 0 

diminishing effect 

both low and high energies 

<5k 7 > 0 
A 7 < 0 

large enhancing effect at low energies 

diminishing effect at high energies, depending on 5 k 7 

<5k 7 < 0 
A 7 < 0 

large enhancing effect at low energies 

small enhancing effect at high energies, depending on <5 k 7 


TABLE III: Disentangling different scenarios based on the energy distribution of the secondary fi. 


muon emerging in the backward region, with with some trace in the high energy - backward region, 
and very low energy forward region. On the other hand, the signature of 5 k 7 will be seen more 
visibly in the forward direction for high energy /i. In the case of low energy fx, the effect of & 7 is 
somewhat independent of the angle, except for the extreme forward region. The reason for these 
complex dependence is a result of the complex interdependence of the kinematic variables and the 
parameters, the origin of which are difficult to decipher,. 

Taking slices from the above surface plots, the energy distributions at fixed angles of 60° and 
150° are presented in Fig. corresponding to center-of-mass energy of 1000 GeV. These angles 
represent the typical behaviour in the low and high angle regions, respectively. Quantifying these 
effects, the number of events in the chosen energy bins for the above angles, corresponding to the 
SM case as well as different combinations of the anomalous couplings, along with the percentage 
deviations is given in Table [TV) The energy bins are chosen to maximize the effect, in each case. 
In the forward region with high energy muon, where the effect of A 7 is small, whereas the 5re 7 has 
large effects has a maximum of 12% in the case of y/s = 500 GeV, which remains more or less 
at the same level (10%) in the case of yfs = 1000 GeV. On the other hand, the presence of A 7 
has much larger impact, with 33% when <5 k 7 is zero at yfs = 500 GeV, which is increased to 41% 
along with positive values of <5 k 7 and decreased to 21% for its negative values. These effects can 
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FIG. 10: The energy-angle double distribution of the secondary muon (along the vertical axis: ^ 

(pb)) at yfs = 500 GeV. Top left is the SM case. Other cases are deviations from the SM; <5 k 7 = 0, A 7 = 
—0.059 (top right); <5 k 7 = +0.069, A 7 = 0 (bottom left); <5 k 7 = —0.072, A 7 = 0 (bottom right). 




y/s = 500 GeV 

y/s = 1000 GeV 



^ = 

60° 

= 

150° 

= 

0 

o 

CD 

= 

150° 



x„ = 0.95 - 1.0 

= 0.056 - 0.106 

Xf, = 0.95 - 1.0 

= 0.056 - 0.106 

Sk~, 

A^-y 

Nevents 

A [%] 

■Nevents 

A [%] 

■Nevents 

A [%] 

Nevents 

A [%] 

SM Case 

13500 

- 

6200 

- 

14750 

0 

6000 

0 

0.066 

0 

14550 

8 

6750 

9 

15900 

8 

7000 

17 

-0.072 

0 

12450 

-8 

5650 

-9 

13600 

-8 

5000 

-17 

0 

0.017 

13750 

2 

5750 

-7 

15000 

2 

4000 

-33 

0 

-0.059 

13000 

-4 

8250 

33 

14500 

-2 

14500 

141 

0.066 

0.017 

14750 

9 

6250 

1 

16225 

10 

4750 

-21 

0.066 

-0.059 

14000 

4 

8750 

41 

15600 

6 

15500 

158 

-0.072 

0.017 

12650 

-6 

5250 

-15 

13850 

-6 

3000 

-50 

-0.072 

-0.059 

11900 

-12 

7500 

21 

13300 

-10 

13400 

123 


TABLE IV: The number of events within specified a+ bins, and for fixed for different combinations of 
<5 k 7 and A 7 at yfs = 500 GeV and 1000 GeV, along with the corresponding percentage deviation (A) from 
the SM case. An integrated luminosity of 100fb _1 is considered. 


be increased to more than 4 times for yfs = 1000 GeV. 


Focusing on the low energy region, we consider variation with 6 M for fixed = 0.2 (as a typical 


value), and present difference scenarios in Fig. 12 While the absence of A 7 and its positive values 
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FIG. 11: The energy distribution of the secondary muon for different combinations of <5 k 7 and A 7 , consid¬ 
ered at muon angles of 6 M = 60° and 150°, at y/s = 1000 GeV. 
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FIG. 12: The energy-angle distribution of the secondary muon for different combinations of 5 k 7 and A 7 , 
considered at fixed muon energy parameter ( x at yfs of 500 GeV and 1000 GeV. Colour coding is the 
same in all plots. 


for different <5 k 7 values have moderate impact, especially in the backward region, the negative 
values of A 7 could produce significantly different distributions. In order to quantify the effects, we 
have computed the forward-backward asymmetry in each of these cases, and tabulated in Table [V] 
The deviations are evidently more pronounced in all cases, and much larger in the A 7 < 0 case. 
Note that owing to large cross-section, even limiting to a small energy bin of = 0.2 around 
the fj, energy considered will lead to about 40000 events, when integrated over the angle. Thus, a 
change in the forward-backward asymmetry to the tune of 200% is certainly detectable. Coming 
to the reach on the value of the couplings achieved through these observations, we note that the 
asymmetry in the SM case is about —19%, and +40% at yfs = 500 GeV and 1000 GeV, respectively. 
We summarize the extend of asymmetry measurement in Table |VI[ considering purely statistical 
uncertainty. This suggests that deviations from the asymmetry of about 5 % could be expected 
at 3a level for . Reading it along side Table [V] indicates a reach of about an order of magnitude 
better than that the value considered there, which means about A 7 ~ —0.003, assuming a linear 
dependence. 
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<5k 7 

Ay 

A Afb{%) 



yfs = 500 GeV 

a/s — 

1000 GeV 

0.066 

0 

-4 

-33 

-0.072 

0 

5 

48 

0 

0.017 

-34 

72 

0 

-0.059 

92 

-205 

0.066 

0.017 

-36 

29 

0.066 

-0.059 

83 

-213 

-0.072 

0.017 

-32 

135 

-0.072 

-0.059 

102 

-194 

SM case; Afb = 

-0.185 

0.398 


TABLE V: Observed forward-backward asymmetry and its deviation from the SM with a fixed muon 
energy corresponding to ay = 0.2, at y/s = 500 and 1000 GeV. 



n f 

N b 

A FB % 

A A fb% 

reach on A 7 
at (3 x A A fb ) level 

SM 

W ± 'JW 

18000 

[17866,18134] 

26000 

[25839,26161] 

18.5 

[17.5,18.8] 

[-5.0,+1.6] 

-0.003 


TABLE VI: Number of events in the forward {N F ) and backward ( N B ) hemispheres, and the corre¬ 
sponding forward-backward asymmetry ( A FB ) in the case of the SM with a fixed muon energy of 50 GeV 
(corresponding to x M = 0.2) considering a bin of Aay = 0.2 (A = 50 GeV), at yfs = 500 GeV and 
luminosity of 100 fb -1 . 


IV. SUMMARY AND CONCLUSIONS 

With the W decaying into muons, tthe single W production in ey collisions has a spectacular 
final state of one single muon, and missing energy. The process is well suited to study the WW 7 
coupling. It may be noted that the process has the advantage that, unlike other relevant process 
like the W pair production at e + e~ collisions, it is devoid of effects from other anomalous couplings 
like WWZ. Considering this process, angular distribution of secondary muons in the lab frame is 
easily constructed, and is readily available in the literature. At the same time, such is not the case 
of its energy distribution. Here we have presented a semi-analytical way to explore the secondary 
lepton energy-angle distributions in ey —» vW with W —> Iv. The advantage of such an observable 
in analyzing the SM case and probing possible new physics effects is demonstrated. Variables being 
defined in the lab frame, are directly used to probe different kinematic regions, so as to explore 
the sensitivity of the anomalous couplings. 

We have derived possible limits on the anomalous couplings Sk 7 and A 7 , which may be obtained 
from the cross-section measurements. It is shown that, assuming the other couplings is absent, 3<r 
level of cross-section can probe Sn-y to the level of ±0.004, and constrain A 7 to —0.013 < A 7 < 
+0.01, for a moderate luminosity of 100 flu 1 at a center-of-mass energy of 500 GeV. The large 
influence of each of these couplings on the sensitivity of the other is quite visible in the derivable 
limits. Scanning the two parameter plane the 3 a bound is presented, showing that i5k 7 could 
be bounded to small regions for particular values of A 7 . The angular distributions along with 
forward-backward asymmetry improve the situation with better sensitivity. With the large cross- 
section, the statistics is very good even for 100 fb -1 integrated luminosity, which could measure 
the forward-backward asymmetry to a few percent level. It is found that this could be utilized to 
improve the limits on A 7 , whereas not so in the case of Sk^. Exploring different kinematic regions, 
and studying suitably constructed observables would be able to distinguish different scenarios 
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involving the couplings. We studied the energy and angle distributions of the n separately and 
together to understand how these could be fruitfully employed for this purpose. The energy-angle 
double distribution enables one to identify the kinematic regions sensitive to different couplings, 
and their signs. We have found that focusing on the low energy fi in the backward direction are 
more sensitive to the case of A 7 < 0, whereas the <5 k 7 does not show distinguishable preference to 
the kinematic regions. Our analysis shows that the case of (i) A = 0, or non-zero and positive, (ii) 
non-zero and negative can be distinguished this way, as explained in Table [TTT| On the other hand, 
the sign of <5re 7 can be distinguished in the absence of A 7 , whereas it is hard to do this if A 7 is 
non-zero. At \/s = 500 GeV, with 100 fb' 1 luminosity, the bound on the A 7 , when it is negative 
could be improved by about a factor of 2 by considering the forward-backward asymmetry at a 
fixed energy (25 M ), compared to the case when the energy is integrated out. 

We are aware that the present considerations of the ILC do not favour a photon collider option 
as priority. At the same time, we hope that this study will help the case for this option at ILC, and 
expect it to be available eventually. Meanwhile, although more accurate estimates incorporating 
collider and detector effects, and beam polarization effects are needed to make more quantitative 
conclusions, we hope our analyses has shown that the leptonic decay of the W produced in ey colli¬ 
sion is one of the neat and simple process to probe the WlVy coupling without the complications of 
other couplings. Further, the process has the potential to probe the couplings to a few per-mil level 
with suitably constructed observables, which could also be used to distinguish different scenarios 
involving the two parameters, Sk 7 and A 7 . 
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Appendix: Photon luminosity distribution 


The colliding photons in a realistic electron-photon collider does not have a fixed energy, rather the 
beam will have distribution of photons with energy varying over an allowed range (which depends 
on the initial electron and laser photon energies among other things). In such colliders, the cross- 
section and other observables should, therefore, be properly folded with a luminosity distribution 
function to get the measurable quantities, as is done in Eq. [12] 

At ILC high energy, high luminosity photon beam is obtained by Compton backscattering of 
low energy, high intensity laser beam off high energy electron beam. Ideal Compton backscattered 
photon spectrum is given by f29Tf8T] 


fy/e{ x ) 

m) 


m 


1 — X + 


1 


-4 


+ 4 


1 — - —— 

e c 2 


l — x £(1 — x) £ 2 (1 — x) 2 _ 
1 8 1 


In (1 + 0 


Z 2 ( 1 + 0 2 


( 11 ) 


where x = -)+, with E e the energy of the initial electron and ui the energy of the scattered photon. 
x thus gives the fraction of the electron energy carried by the scattered photon. Dependence of 
the distributions on the initial laser photon energy (wo) comes through £ ss , where m e is the 

electron mass. The maximum value of x is x max = It is, but not possible to increase ojq and 
E e to any value to get larger x max . It is found that for £ beyond ~ 4.8, conversion efficiency drops 
down drastically due to e + e _ pair production between the laser photons and the backscattered 
photons, setting an absolute upper limit on i w 0.83. This value essentially means that with an 
electron beam of energy E e = 250 GeV, we can effectively go up to ujq ~ 1.26 eV . 

For an ey collider, the luminosity factor (Eq. 12) is the same as the above energy spectrum: 


d/l 7 / e (x) 

dx 




( 12 ) 


In a realistic collider, one need to also worry about many detailed aspects, like the non-linear 


effects making the actual photon spectrum deviating from the ideal case in Eq. 11 the polarization 
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of the hard photon. The polarization (P 7 ) of the scattered photon itself depend on the initial 
electron (P e ) and the laser beam polarization (Pi). In the high energy region, it is noted that P 7 = 
— P e [29U3T] , This could be used to the advantage of the physics studies being considered. In our 
analysis, we have not gone into these details, we have assumed unpolarized initial electron and laser 
beams, which produces an unpolarized high energy photon beam (to a very good approximation). 

For more details on photon collider one may refer to Ref. PTH - DTH [53j and references therein. 
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